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The (2 + l)-d U(l) quantum link model is a gauge theory, amenable to quantum simulation, with 
a spontaneously broken SO(2) symmetry emerging at a quantum phase transition. Its low-energy 
physics is described by a (2 + l)-d RP(1) effective field theory, perturbed by a dangerously irrelevant 
SO (2) breaking operator, which prevents the interpretation of the emergent pseudo-Goldstone boson 
as a dual photon. At the quantum phase transition, the model mimics some features of deconfined 
quantum criticality, but remains linearly confining. Deconfinement only sets in at high temperature. 



Quantum link models (QLMs) are lattice gauge the- 
ories formulated in terms of discrete quantum degrees 
of freedom. ^(1) and SU(2) QLMs were first con- 
structed by Horn in 1981 pQ, and further investigated in 
[2]. In [3] QLMs were introduced as an alternative non- 
perturbative regularization of Abelian and non-Abelian 
gauge theories, in which ordinary gauge fields emerge 
dynamically from the dimensional reduction of discrete 
quantum link variables. Dimensional reduction of dis- 
crete variables is a generic phenomenon in asymptotically 
free theories, which gives rise to the D-theory formula- 
tion of quantum field theory [4 j. In the D-theory for- 
mulation of 4-d Quantum Chromodynamics (QCD), the 
confining gluon field emerges by dimensional reduction 
from a deconfined Coulomb phase of a (4 + l)-d SU(3) 
QLM [5 j. Chiral quarks arise naturally as domain wall 
fermions located at the two 4-d sides of a (4 + l)-d slab. 
The (2 + l)-d U(l) QLM has also been investigated in 
the context of quantum spin liquids [6|. With staggered 
background charges ±1, it is equivalent to a quantum 
dimer model [7-9 . Furthermore, Kitaev's toric code [10] 
is a Z(2) QLM. In contrast to Wilson's lattice gauge the- 
ory [11], QLMs have a finite-dimensional Hilbert space 
per link, which makes them ideally suited for the con- 
struction of atomic quantum simulators for dynamical 
Abelian [T2UT6] and non-Abelian gauge theories [T7H20] . 
A long-term goal of this research is to quantum simulate 
QCD in the D-theory formulation with ultracold matter, 
in order to address the real-time evolution of strongly in- 
teracting systems in nuclear and particle physics, as well 
as their dynamics at non-zero baryon density. 

In this paper, we investigate the (2+l)-d U (1) QLM, in 
order to demonstrate that, despite its structural simplic- 
ity, it displays highly non- trivial dynamics, and thus is 
ideally suited to demonstrate the power of gauge theory 
quantum simulators. We consider the model with a pla- 
quette coupling J and a Rokhsar-Kivelson (RK) coupling 
A. The phase diagram is sketched in Fig.l. At zero tem- 
perature, the model is confining for A < 1. At finite tem- 
perature T, it has a deconfinement phase transition above 
which there is a massless mode transforming non-trivially 
under the U(l) center symmetry. Due to the Mermin- 
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FIG. 1. [Color online] Schematic sketch of the X-T phase 
diagram. The insets indicate the location of the peaks in the 
probability distribution of the order parameter p(Ma, Mb)- 



Wagner theorem, this is not associated with spontaneous 
symmetry breaking. The deconfinement phase transition 
reaches zero temperature at the RK point, A = 1. At A c 
there is a quantum phase transition which separates two 
phases with spontaneously broken translation symmetry 
[21]. The phase at A < A c has, in addition, a sponta- 
neously broken charge conjugation symmetry. The two 
phases are similar to the columnar and plaquette ordered 
valence bond solid phases in a quantum dimer model [8] , 
which may be separated by a first order phase transition 
[22] or by an intermediate phase [9]. 

As we will see, at A c a spontaneously broken approxi- 
mate global SO (2) symmetry emerges dynamically, giv- 
ing rise to a light pseudo-Goldstone boson. The interface 
that separates the two broken phases on either side of the 
transition manifests itself as a string with fractional elec- 
tric flux |. This raises the question whether the phase 
transition might be a deconfined quantum critical point 
in the sense of [23] [24] , corresponding to a conformal field 
theory with an emergent massless photon and deconfined 
electric charges. Deconfined quantum criticality has first 
been investigated numerically in the J-Q quantum spin 
model [25-27 , between an antiferromagnetic and a va- 
lence bond solid phase, and is still discussed controver- 
sially [28] [29]. It has also been studied in the J-Q [30] 
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and in the J1-J2 model [31-33 on the honeycomb lattice. 
At a deconfined quantum critical point, the instanton- 
like monopole events that cause permanent confinement 
in a (2 + l)-d compact U(l) gauge theory [34j [35] are 
eliminated, because a Z(4)-invariant term that explicitly 
breaks the emergent SO (2) symmetry of the effective ac- 
tion becomes irrelevant. We will see that this is not what 
happens in the (2 + l)-d 17(1) QLM, where the Z(4)- 
invariant term can be tuned to zero. Still, the emergent 
SO (2) symmetry remains weakly explicitly broken by a 
"dangerously irrelevant" operator [36, 37 . This prevents 
the interpretation of the Goldstone boson as an emer- 
gent dual photon. It is more appropriate to think of it 
as an accidentally light Abelian "glueball" . The danger- 
ously irrelevant operator also contributes to the string 
tension and implies that the theory remains confining at 
the phase transition. 

The Hamiltonian of the (2 + l)-d 17(1) QLM is 



H 



-J 



□ 



\(U U + U ] u f 



(1) 



Here U\j = U wx U xy U\ y U} )JZ is a plaquette opera- 
tor formed by quantum links U xy connecting nearest- 
neighbor sites x and y on a 2-d square lattice. A J7(l) 
quantum link U xy = S xy is a raising operator of electric 

flux E = S xy1 constructed from a quantum spin S xy as- 
sociated with the link xy. In Wilson's lattice gauge the- 
ory, where the link variables are classical parallel trans- 
porters, U xy = ex.-p(i(p xy ) G £7(1), taking values in the 
gauge group, and E xy = —id^ xy , the single-link Hilbert 
space is infinite-dimensional. In the U(l) QLM, on the 
other hand, it is just given by a finite-dimensional rep- 
resentation of the embedding algebra SU(2). When one 
chooses spin | on each link, the link Hilbert space is 
just 2-dimensional. The first term in the Hamiltonian 
flips a loop of electric flux, winding around an elemen- 
tary plaquette, and annihilates non-flippable plaquette 
states, while the RK term, proportional to A, counts 
flippable plaquettes. The Hamiltonian commutes with 
the generators, G x = Y,i( E x ,x+i E x-i,x)' of infinitesi- 
mal U(l) gauge transformations. Here i is a unit- vector 
pointing in the i-direction. Physical states are gauge 
invariant, i.e. they obey the Gauss law G x \^/} = 0. 
Besides the gauge symmetry, the QLM also has sev- 
eral global symmetries, including lattice translation in- 
var iance and charge conjugation. Translation invariance 
characterizes each energy eigenstate by its lattice mo- 
mentum p = (pi,P2) £ ] — K^Tr] 2 . Charge conjugation 
replaces U xy by U% y and E xy by —E xy , and charac- 
terizes each eigenstate by its charge conjugation parity 
C = =b. Another important global symmetry is the cen- 
ter symmetry associated with "large" gauge transforma- 
tions. The £7(1) QLM defined on a periodic volume has 
super-selection sectors characterized by wrapping electric 
fluxes that take values in Z/2. On an Li x L2 lattice with 
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FIG. 2. [Color online] a) Energy gaps of the lowest states on 
the 6x6 lattice as a function of X. b) Global fit of L1L2- 
dependence of X pc , \' p t , and X' pc , that yields X c = —0.359(5). 



periodic boundary conditions, the generators of the U(l) 
center symmetry are Ei = j- ^2 X E x x+ ^. They commute 
with the Hamiltonian, but cannot be expressed through 
"small" periodic gauge transformations G x . 

We have performed exact diagonalization studies of the 
(2 + l)-d 17(1) QLM with S = \ on 4 x 4, 4 x 6, and 
6x6 lattices. The energies of the lowest states are il- 
lustrated in Fig. 2a. For A < 1, the ground state has 
momentum (0,0) and is even under charge conjugation 
(i.e. C = +). For A < A c the first excited state has 
quantum numbers C = — , p = (tt, it). Its energy gap 
to the ground state, E_ ~ exp(— <7_LiL 2 ), decreases ex- 
ponentially with the volume L1L2, thus indicating the 
spontaneous breakdown of charge conjugation C and the 
translation T by one lattice spacing (in either direction). 
For A > A c , another state \C = +,p = (7r,7r)) degen- 
erates with the ground state, i.e. E+ ~ exp(— cr + LiL 2 ), 
indicating that C is now restored, while T remains spon- 
taneously broken. The crossing of the two excited energy 
levels, E_(X pc ) = £? + (A pc ), defines a volume-dependent 
pseudo-critical coupling X pc . The next zero-momentum 
excited states, \C = ±,p = (0,0)) with energies E±, 
cross twice near the critical point at two pseudo-critical 
couplings X p + and X p ~, i.e. E'_(X' p f) = E+(\' p f). As il- 
lustrated in Fig. 2b, the couplings A pc , A^, and X p ~ all 
approach A c = —0.359(5) in the infinite volume limit. 

The different symmetry breaking patterns are distin- 
guished by two order parameters, Ma and M#, associ- 
ated with the even and odd dual sublattices A and B. 
A configuration of quantum height variables h~ = 0,1, 
h~ = d=|, located at the dual sites x = (x\ + ^,£2 + |), 
is associated with a flux configuration E x = [h~ — 

fc^_ i _ 2 ]mod2 = ±|, X,X' e {A, B}. The two or- 
der parameters are given by Mx = ^2 x ex s ^^^ where 

S§ = (_l)(^i-^)/2 and S B = (_ 1 yx 1 -x 2 + l)/2 m Under C 

and T they transform as C M A M A , C M B = -M B , 
T M A -M B , t M b = M A - It should be pointed out 
that ±(Ma, Mb) represent the same physical configura- 
tion, because shifting the height variables to h~ (t)' = 
[h~(t) + l]mod2 leaves the electric flux configuration un- 
changed. We have performed quantum Monte Carlo sim- 
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FIG. 3. [Color online] Probability distribution p(M a, Mb) for 
A = — 1, A C; and at T = (a,b,c), and A = 0, T > T c (d). 



illations with an efficient newly developed cluster algo- 
rithm, that will be described elsewhere. The algorithm 
has been used to determine the probability distribution 
p(Ma,Mb) of the two order parameters Ma and Mb 
shown in Fig. 3 at A = — 1, A c , and for L\ = L2 = 24a, 
which reveals an emergent spontaneously broken SO (2) 
symmetry at the quantum phase transition. 

The low-energy effective theory describing the vicinity 
of the quantum phase transition is formulated in terms 
of a unit- vector field e(x) — (cos ip(x), sin <p(x)) repre- 
senting the direction of (Ma, Mb)- Since (Ma, Mb) and 
— (Ma, Mb) are indistinguishable, the effective theory is 
a (2 + l)-d RP(1) model. Thus, only those states that 
are invariant against a sign-change of e(x) belong to the 
physical Hilbert space. Introducing ds = d ctl the corre- 
sponding Euclidean effective action is 



S[<p]: 



d 6 x 



d fI (pd fI (p+5cos 2 (2(p)+£cos 4 (2(p) . (2) 



2~ M 



Here p is the spin stiffness and c is the velocity of an 
emergent pseudo-Goldstone boson. 5 + e measures the 
deviation from the phase transition. The 5-term explic- 
itly breaks the emergent SO (2) symmetry to a Z(4) sub- 
group and gives rise to a small Goldstone boson mass 
Mc — 2yj2\b\j p. Even when the relevant ^-term is tuned 
to zero, the dangerously irrelevant £-term still explic- 
itly breaks the SO (2) symmetry. It is natural to de- 
fine the dual field F fjLU (x) = \e pv p d p tp{x) . Since ip(x) is 
well-defined only up to multiples of 7r, vortices and half- 
vortices in the order parameter field manifest themselves 
as charges. The electric charge contained in a spatial 
region Q is given by twice the vortex number 



Qn = / d 2 x diF 0i = - [ dai Sijdjcp e 



(3) 
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Note that a charge 1 corresponds to a half-vortex, which 



FIG. 4. [Color online] Phase diagram as a function of 5 and 
e. The insets indicate the location of the peaks in the dis- 
tribution p(Ma, Mb)- The fat and dashed lines are first and 
second order phase transitions, respectively. The curved line 
indicates a possible path taken in the QLM when varying A. 



is allowed because e(x) and —e(x) are physically equiva- 
lent. While the flux of F pv correctly represents the con- 
served charges of the U(l) center symmetry, F pv should 
not be mistaken for a dual massless photon. This in- 
terpretation would require an exact 50(2) symmetry, 
at least in the infrared. Due to the £-term and other 
higher order symmetry breaking terms, this would re- 
quire a large amount of fine-tuning. 

By applying the Ginsburg-Landau- Wilson paradigm to 
the 5- and £-terms, in mean field theory one obtains the 
phase diagram of Fig. 4. The two phases realized in the 
QLM both have four peaks in the order parameter dis- 
tribution p(Ma, Mb), and are separated by a weak first 
order phase transition. In addition, there is an interme- 
diate phase with eight peaks (whose analog may be re- 
alized in the quantum dimer model [9 ), separated from 
the other phases by second order phase transitions [36] . 
If one would fine-tune to these transitions, the Goldstone 
boson would become exactly massless. Even then it could 
not be interpreted as a dual photon, because the 50(2) 
symmetry would still remain explicitly broken. 

The effective theory predicts a finite-volume rotor 
spectrum (E m = m 2 c 2 /(2pL\L,2) when 5 = e = 0) with 
m = 0, ±2, ±4, . . . States with odd values of m are ex- 
cluded because they are not invariant against a sign- 
change of e(x). The quantum numbers of the states 
with m = 0,±2,±4 correspond to C = +,p = (0,0), 
C = ±,_p = (7r,7r), C = ±,p = (0,0), respectively. The 
effective theory is in quantitative agreement with the ex- 
act diagonalization study (cf. Fig. 2b). A global fit of the 
energy spectrum yields A c = —0.359(5), p = 0.45(3) J, 
c = 1.5(l)Ja, 5 C = —e c = 0.01(l)J/a 2 . A more pre- 
cise determination of the low-energy constants and of A c , 
based on high accuracy Monte Carlo simulations using 
the cluster algorithm, will be presented elsewhere. 

Away from the critical point, the 5-term gives rise to 
two distinct coexisting phases that are related by C and T 
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FIG. 5. [Color online] a) Interface at A — >> — oo. The fat line 
is a flux \ string, b) Potential between two static charges ±2 
separated by the distance along a lattice diagonal, for 

A -1, A C; and 0, at T = 0, and at A = A c /or T 2 J. 



for A < A c , and by T for A > A c . As illustrated schemat- 
ically in Fig. 5a for A — » — oo, the interface that sepa- 
rates the two phases represents a string of electric flux ^ . 
Its interface tension, = \/2|5|p (for e = 0), which 
plays the role of a string tension, would vanish at the 
phase transition if there was no £-term. While the e- 
term would simply be irrelevant at a critical point, here 
it is dangerously irrelevant. Taking it into account, the 
string tension never vanishes, and is always of order yfpe. 
Indeed, the potential (shifted by a constant) between two 
static charges ±2, illustrated in Fig. 5b, shows linear con- 
finement at large distances, even at the phase transition, 
albeit with a small string tension cr 2 = 0.201(2) J/ a (com- 
pared to (72 = 1.97(1) J/ a at A = —1). This shows explic- 
itly that the phase transition is not a deconfined quan- 
tum critical point. The energy density —J{Uu + U^) in 
the presence of two charges ±2 is illustrated in Fig. 6a- 
d. The flux string connecting the charges separates into 
four strands of flux \ that repel each other. In accordance 
with the effective theory, the interior of the strands con- 
sists of the phase that is stable on the other side of the 
transition. Near A c the flux string undergoes topology 
change by wrapping one strand over the periodic bound- 
ary and materializing an additional strand at the edge 
of the system, whose interior then expands to become 
the new bulk phase (cf. Fig. 6b). Viewed as interfaces 
separating bulk phases, the strands display the universal 
phenomenon of complete wetting. 

Finally, we have studied the system at finite temper- 
ature. The electric flux susceptibility, (Ef) is non-zero 
at T > T c , indicating a massless mode that transforms 
non-trivially under the U(l) center symmetry, giving 
rise to a logarithmic charge- anti-charge potential (cf. 
Fig. 5b). Hence, the "deconfined" phase no longer has lin- 
ear, but still has logarithmic confinement. As illustrated 
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FIG. 6. [Color online] Energy density —J(U\j + £/"□) in the 
presence of two charges ±2 for A = — 1 (a), A = A c (b), and 
A = (c) atT = 0, as well as for A = at T > T c (d). 



in Fig.6d, the flux then spreads out and no longer forms 
a string. Interestingly, the shift symmetry T remains 
spontaneously broken at high temperature (cf. Fig. 3d). 
Actually, in the deconfined phase yet another SO (2) sym- 
metry emerges, which originates from the Gauss law. 

In conclusion, we have observed an emergent SO (2) 
symmetry with an associated pseudo-Goldstone boson in 
the (2+l)-d U(l) QLM. Interfaces separating phases with 
spontaneously broken C or T symmetry manifest them- 
selves as strings carrying fractional electric flux |. Al- 
though the model displays certain features of deconfined 
quantum critical points, a dangerously irrelevant opera- 
tor leads to small explicit SO (2) breaking. This prevents 
the interpretation of the emergent Goldstone boson as 
a massless photon, and implies a non-zero string ten- 
sion also at the phase transition. It remains to be seen 
whether phenomena, similar to the ones observed in the 
QLM, may masquerade as deconfined quantum critical- 
ity in other models as well. Once the (2 + l)-d U(l) 
QLM is realized in ultracold matter experiments, its rich 
dynamics will become accessible to quantum simulation. 
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